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SECTION  I 

ASYMPTOTIC  FORMULAE  FOR  THE  DISTRIBUTION  OF  HOTELLING'S  TRACE 
FOR  TESTS  OF  EQUALITY  OF  TWO  COVARIANCE  MATRICES 

1.  INTRODUCTION 

Let  m Sj  and  n S2  be  independently  distributed  W(m,p,Zj^)  and  W(n,p,Z2)» 

respectively.  Chattopadhyay  and  Pillai  [ 1]  have  given  asymptotic  expansions 

for  the  c.d.f.  and  percentiles  of  T = in  TrS^S^^  up  to  terms  of  order  n'^  in 

which  the  noncentrality  was  denoted  by  (F)  = tr  F = tr(B  ^A  - I) /'the 

deviation  matrix",  where  B = and  A = In  their  paper^  terms  involving 

fijfki* where  f^j  is  the  (i,j)  element  of  F,  have  been  neglected.  These  terms 

are  taken  into  consideration  in  the  section,  and  noncentrality  is  expressed  in 

the  form  (Fs)  = tr  F^.  Table  I gives  tabulations  to  show  the  importance  of 

these  terms.  Furthermore,  Chattopadhyay-Pillai  (denoted  by  C-P 

2 

hereafter)  expansions  are  extended  to  terms  of  order  1/n  . It  may  be  noted 

here  that  T = nU^**^,  where  U^**^  is  the  statistic  studied  by  Pillai  [2]  for 

the  test  of  « Z2,  and  the  power  of  this  test  against  alternatives  of  a 

one-sided  nature  was  studied  by  Pillai  and  Jayachandran  [3]  for  p>2. 

Recently  the  exact  non-null  distribution  of  Hotelling's  trace  and  tabulations 

of  the  power  of  the  same  test  for  small  and  large  deviationlof  the  parameters 

were  studied  by  Pillai  and  Sudjana  [4]  for  p=3  and  m>4.  Some  power 

_2 

tabulations  are  presented  in  Table  I up  to  terms  of  order  n , which  show 
extremely  good  accuracy  compared  to  the  exact  values  given  by  Pillai  and 
Sudjana.  Some  additional  tabulations  of  powers  are  also  presented  for  p«4. 

2.  THE  METHOD  OF  ASYMPTOTIC  EXPANSION 

The  notations  here  as  well  as  in  the  rest  of  the  section  follow 
those  of  [1  ] and  [5].  In  order  to  describe  the  method  we  will  first  derive 
an  asymptotic  expansion  for  the  percentiles  of  T using  which  we  will  further 


1 


I 


obtain  that  of  the  c.d.  f.  of  T.  It  is  well  known  [6  ] that  the  statistic 

1ft  2 2 

y ■ ■ Tr  Sj  A can  ba  writtan  as  £'  X^Xj  (■)*  whara  Xj(a)  '•  indapandent 
central  chi-square  variables  with  a d.f.  and  X^'s,  j ■ l,...,p,  are  the 
characteristic  roots  of  U ■ A 

m m 

Let  G(0)  ■ Pr{a  Tr  S^A  < 26). 

Now  note  that 

Pr{a  Tr  SjB  < 20)  - G^(0)  - [r(p)]'VJ  iV’^dt,  (1) 

where  p ■ ap/2.  In  G(6),  as  a first  approxiaation,  for  large  n we  aay 
replace  A~^  by  S,  and  consider 

G(0)  - Pr{a  Tr  Sj  S"^  < 20).  (2) 

Now  we  nay.  as  suggested  in  [ 5 ].  obtain  a function  h(S2)  in  the  elenents  of 
such  that 

G(b)  - Pr{a  Tr  < 2h(S2)), 

and  then  write  h(S2)  as  a series  with  the  first  tera  being  a linear  function 
of  chi-square  variables  and  successive  texas  of  decreasing  order  of  nagnitude. 

We  get 


Pr{a  Tr  SjS“^  < 2h(Sj))  - / Pr{a  Tr  SjS’^  < 2h(S2)IS2)  Pr{dS2),  (3) 

where  Pr{dS2)  is  the  probability  elenent  of  the  central  Nishart  distribution 
of  $2  and  R is  the  donain  of  integration  of  $2.  Now  we  nay  expand 
Pr{®  Tr  ^^^2  — 2h(S2)|S2)  about  an  origin  (®jj»®22’''’**^pp*^12’**’*^p-l  p^ 
in  a Taylor  series,  idiere 


.-1 


(o^j)  i.  J ■ l....,p. 


(4) 


Thus 


Pr{n  Tr  SjS'^  < 2h (82)182) 


■t«xp[  I (8jj  - Ojj)  ])  Pr  {a  Tr  8^A  ^ 2h(A“^)) 
1<J-1 
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• M 


{exp[Tr(S,  - A‘^)3])  Pr  {■  Tr  S.A  < 2h(A"^)), 

m «X«  m 


where 


3 

■5o]J 

3 

1 3 

1 3 

L 

^ *«12  

3 

1 3 

V 

*"21 
1 3 

*“22 

1 3 

•••  * *“2p 
3 

* • »"pl 

2 ‘ »‘='p2  *** 

PP 

where  6^^  is  the  Kronecker  delta.  Hence  substitution  of  Eq  (5)  into  Eq  (3) 
and  tern  by  term  integration  for  sufficiently  large  n gives 

G(0)  - / exp[Tr(S,  - A’^)3]  Pr  (a  Tr  S.A  < 2h(A"^)>  Pr{dS,} 


where 


- ® Pr{a  Tr  S.A  < 2h(A‘b>, 

«X  ^ 


® » exp[-Tr  A“^3](r  (n))"^  |a|"^^. 

• P M 


/ exp[Tr(S,  a-  (n/2)AS,]dS 


#w2j  ^*^2 


exp(-Tr  A"^3]|l  - (2/n)  A"^ 


Now  using  [ 7],  we  get 

0 ■ 1 ♦ - 3^^  3^^  ♦ -j 

>.t v«  * I ZVtu“vwV.tV«V’ 

fdiere  £ denotes  the  sunmetlon  over  all  suffixes  each  of  which  ranges 

from  1 to  p.  FVirther,  we  represent  h(S.)  as 


3 


(8) 


h(S^)  - e+  hj(S^)  h2(S^H. 


where  ^^(82)  is  of  order  n"®.  Then  Eq  (8)  may  be  substituted  into 

Prfm  Tr  Taylor's  expansion  we  have 


PrCm  Tr  S A < 2h(A"^)) 


:xi 


» exp[{h-(A"^)  + h,(A"^)+ )d]  Pr(m  Tr  S.A  < 20) 

- [1  + {h.  (A‘^)  + h,(A"b+....)D  + -JCh-CA"^)  + 


-U 


^{*  Tr  ^A  < 20), 


where 
Hence  we  get 


„ a 

““■Se 


G(0)  ■ Cl  + ” St  <t  +”«»{^2oo<t  d & d + 

n re  tu  St  ur  ^2  h rs  tu  vw  st  uv  wr  ^ 

i ^ ®rs®tu°v«Vtt®urVyv’  + “<«''>  3tl  + 1i,(*'Sd  + 

{h,(*‘‘)D  + 4f(A*‘)D*J  + 0(n'’)]  Pr{«  Tr  S,A  < 29). 


Now, equating  terms  of  successive  order  [ 1]>  we  have 


{hi (A  ^)D  + ^ L Pr  (.  Tr  S^A  < 20)  - 0, 

Ch2(A"^)D  + V(A’bD* 

■"n^Wtu  »■*)!>+  2h{‘*>  (A-*)  9 D 

^ * SM  U* 


+ d D)  +S^z  a a c d & d 

‘ ~ St  ur  2 rs  tu  vw  st  nv  i 


3n 


wr 


^ ^ ^ VtuVxy^t^rVyv^  20)  - 0, 


(9) 

(10) 


(11) 


(12) 


isittUiiiiimiiiiiiiieieiid^^ 


\ 

J 


nd  ..  <*,  h<“>(A-‘)  . S„h,(A*>)  tnd  h<“-”>(A-‘)  . a 9 h (A**) 

Hence  to  evaluete  h (a‘^)  and  h,(A‘^)  «•  have  to  find 

^•t^ur  ^ ^t^uv^wr  ^ ^ 20), . . . • tor  this  purpose  we 

use  perturbation  technique  [ 8 ] . Let 

J - Pr{m  Tr  ^(a"^  + «)"^  < 20), 


Idler.  «(p»p)  t,  , ,jBi«tric  aetriz  .ufftcleetly  do.,  to  OCpzp).  ».  I.vlot'. 

M ^ 

theorem  we  get 


J - {I  * E «r.»„  * ir  ^ ‘r.*t«*r.*tu  * fi  = Vr.Vw  * 


41  ^ *M*tu*vu^*y^Bi^tu^vir^xy*" ^ Pr{m  It  ^A  < 20). 


C13) 


Also  by  definition  we  get 

|b|“'2 

^ " ; „ (mn)/?  f ® Y ^'3  df, 

(2n)'"P^'2  .'n  2 ...» 

where  « • W',  Y(paO  and  R;{Y:»  Tr  S^(a’K  < 29).  Now  let  r(p«p) 

be  a nonalngular  matrix  auch  that 

T r'  B r - I(pKp)  - D-, 

A M ^ M M ^1 

and 

2 V £•  £<P*P). 

for  e(pxp)  sufficiently  cleee  to  O(pxp)  and  D«  - dUg  (7), T)  ) . This 

.1  * P 

Is  possible  as  B and  A ^ are  p.d. 

Let 

Y(pxm)  ■ r(pxp)  Z(pxm}. 

Mm 


5 


Thm 


Cp(e>. 


where  p - ep/2  and  E la  an  operator  such  that  E 6 (0)  - G (6) . Row  let 

P P+1 

E - A 1. 

Then 

|l  - D e|/|i  - D J - |l  - - MI/|i  - D I 

• |l  - [b‘^(A‘^  + «)“^  - I]A| 

M ^ aw  ^ 

- |l  - XA|,  (aay). 

Hence 

J - |l  - G (0) 

MM  p 

- exp[(-m|2)  log|l  - XA|]  G,(0). 

M ^ P 

Now,lf  b'^A  - I + P such  that|ch.(P)|  <1,  1 - l,...,p,  then  for  «(p*p) 
aufflclently  dote  to  O(pxp)  we  get  |ch.(X)|  <1,  I • l,...,p,  and 

J - exp{f  Tr  XA  + a Tr  X^A*  + f Tr  X^A^  + f Tr  X V + ) G (9) 

£ ^ ^ ^ Oav  P 


- Cl  + f Tr  XA  + {S  Tr  X^  + f-  (Tr  X)*)  A^ 

2 3 

+ {f  Tr  X^  + f-  (Tr  X)(Tr  X^)  + fr  (Tr  X)^)  A^ 

2 2 2 
+ {f  Tr  X^  + ^(Tr  X)  (Tr  X^)  + fr  (Tr  X^) 

V av  ^ a^  3Z  av 


+ ^ (Tr  X)^(Tr  X^  + f«7  (Tr  X)V+ ] G (0). 


32 


384 


(14) 


Now, using  Taylor's  expansion  for  A ^ <f  c,  we  can  represent  X by 

^ awa  - 

X - b"^(a“^  + e)"^-  I 

M M M ^ ^ 

- + I ^ 

^ ^ rPa^P  aM  #; 


- b"^(I  + Z e a A - I 

M #jc*  ^ ** 

- b"^(I  - Z (A  A)  + Z e « (a;;^  AHA^^A) 

- Z c c,  t (A“^  A)  (A"^  A)  (A“^  A) 

r«  tu  w or«  ^ JCu  ^ ^ 

+ Z c c,  e (A‘^  A)  (A"^  A)  (A*^  A,  (a"^  A) 
rs  tu  VH  xy  ^T8  M ^ ^vw  J J 


.)  A - I 


- (b"^A-I)  - Z e (b"^A)(A“1  A) 

+ I s (b"^A)(A“^  A)(A'^  A) 
rs  tu  ^ \rs  J \jn  ^ 

•Sc  c^  «.^<b"^A)  (A'^  A)  (a“^  A)  (A*^  A) 


+ ^ 

^ M ^ ^ ^xy  M 


(15) 


where  A*^  Is  pxp  oatrlx  obtained  by  operating  d on  A*^;  l.e..  It  has  Its 
^s  rs  M 

(l,J)-JJll  element  as  ^ ^***  notations 


Tr(A'^  A)  - (rs), 

xJCB  ~ 

Tr(Ali  A)(a"^  A)  - (rsitu), 

ama»9  a^U  aM 

Tr(A“i  A)(a“^  AHA^^A)  • (rsitulw), 

a^V  av  aJCU  aM  a^vw  av 

Tr(^i  AHA’i  A)(a;2  ^>  " (r»Uu|vw|xy) . 

aws&B  ^ ^ av^W  av  a^y  av 

Tr  F (A’J  A)(a“^  A)  - (f| rsitu), 

^ amXv  Aa  a^U  aM  av 

Tr(B’^  A)(A"i  AHaIJ  A)  - (I  + pjrsltu). 

av  M awA  B aha  aJCU  ah»  aM  aha 

Tr(B"^  A)(A*i  A)(a’^  A)(a;2  A)  - (I  + F|rs|tulvw)  , 

^ M ^ a^U  SM  ah^W  a^  av  av 

Tr  F - (F) 

av  aM 

Tr  F^  ■ (f|  F)  . . . or  alternatively, 

aw  /w  aw 

Tr  F*  - (Fs) . . .etc. 

aw  aw 

7 


- ■ - - 


and  •ubstitutlng  Eq  (15)  into  Eq  (14) , we  get  from  term  by  term 
comparison  between  two  expansions  of  J,  Eqs  (13)  and  (14)  after  substitution, 
the  following: 


I 


j 

I 

1 

1 


d Pr(n  Tr  S.A  < 29)  - -{(f)  (I+F|r«)A  + [(f)  (F|l+r|r«) 

2 

+ (?-)  (F)  (I+f|  re)  + [ (f)  (f|  F j I+f|  re) 

2 3 

+ (f-)  {(F|F)(I+F|re)  + (2)(F)(F|l+F|re))  + ^(F)^I+F|rs)]A3 

2 

+ r<f)(FlFjFll+Flre)  + (fr)  {(6)  (F)  (F|F|l+F|re) 

i-  m « M M M **»  ~ M M M M 

+ (2)(P|F|F)(1+Flre)  + (3)  (F|F)(P|l+F|ra)) 

+ (^)C(F)^(F|l+F|re)  + (F|F)(F)(I+F|rs))  + (^)  (F)^(I+F|re)]d^ 

10  0^  0^  0^  0^  ^ 0^  0m  0^  0^  •50"?  ^ ^ 0m 

^ } Cp(9), 

SleilUrly  expreeelone  for  Pr(n  Tr  S^A  < 20),  Pr(«  Tr  S^A  < 29) 

^re^tu^w^xy  ^ JaI-  “•  available  in  [9  ]. 


3.  AN  ASYMPTOTIC  EXPANSION  FOR  PERCENTILES 
OF  T-m  Tr  S^S'^  UP  TO  0(n'^). 

Now  recall  that  AGp(6)  = -Egp(0),  where  gp(6)  « [r(p)]"^  e'® 
and  p » mp/2.  Also  we  note  here  that 


e‘  i^(e)  - 


9* 

p(p+-l)  • • • (pfl-l) 


*p<®> 


(16) 


Thue.  It  la  poeelble  to  write  PrCei  Tr  S^A  < 20]  in  the  following  form: 


*.t*ur  fliS 


4 

£ 


(17) 


vhere  J - l 4|«re  also  available  in  [9].  Also,  we  note 


^ Vtu  • 9 P<P  + 1). 


t.s.t.u 


\ Vtu  <»t)(ur)  - p.  E a (at)  - p, 


t #u 


“ " ^ , (P|«t)(Br)  - (P), 

r,s,t,u  ” ~ 

^ ^rs^^tu  " (F)(p+l)/2, 

r.s.t.u  ~ ~ 


I • t • • 


A.  . ch.ek  f«r  the  ebm  reUtlonehlpe,  let  F(|a,)  - I(p»p),  ihu.  V ehould 

•hF  mt 

be  e.p«t  to  p,id.teh  U the  velue  of  E a„<t^(et)  (ur) . StaUerly  V end  w 

•hould  eqtiel  | p(prt),»blcb  le  the  velue  of  E <7^e  (,t|ur) . With 

the  eld  of  these  ceeulte  » eveluete  Aj'e,  J - i 

over  all  subscripts  s,  t,  u,  r. 

Hou  bj-  uelus  E,s  (11),  (16),  (17)  and  (18),  and  putting  2e.y  we  get 


Cc'(e)]->, 


where 


^ ‘^rs'^tu  N ^ ^"P>  <"P^2)  ....  (mpf2 J-2) , j 

*»Stt»U  ■* 

have  been  evaluated  and  are  available  in  [ 9], 


^ find  Aj  coefficients 


9 


Ih«n  «•  gtt  froB  Eq  (8)  the  following: 


I - ■ Ir  SjS^*  . , + * o,„-S) 

■ i Jj  «/9)tc'(9)]-l  + o(„-^. 


Hence  we  have  the  following  theorea: 

Iheom  1 Ut  ■ Sj  and  n ^ be  independently  diatributed  W(b,p,b'^), 
»(n.p.A->).  r«p.ctlT.ly.  ,nd  i„  |a^  („|  < . , 

H»n  « ..yptotle  «.„n.lon  top  th.  p.rc«.tu.  .f  I . , j,  S,s:‘l.~*l^  Ty" 

Eq  (20). The  following  are  special  cases  of  Eq  20. 


SMta.  «b«  ««  ^ 

(ij)  of  F,  ton.  Ilk.  (£)*.  (K),  („(„) 

Con..,««.tly  tu.p^r  A,.  A^  .«■  T,  ^ll  k.  ..  .k, 

following, 

A,  • (1/4)  ((pfl)  - ■((P)(pfi)/2  ♦ l)  4 
Aj  - (1/4)  ((pfl)  ^ „(i  . 2{f)/p)  • w^P))/(^2)  end 
Aj  - (1/4)  ((2)  (F)  (pfl) /p  + wt(F)(p+l)/2  4 2(F) /p) 

+ (b^/2)  (F)  ) / («p+2)  . 

and  froB  Eq  (19)  ••  get 

v'] 

which  agrees  with  C-P[l  ] to  the  indicated  order  after  siaplification. 
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2.  As  defined  earlier. 


P 2 

y • S X.X.(n) 
j-1  J J 


2 

where  X.(n)'s  ere  Independent  eencrel  ehl-squere  variables  with  ■ d.f.  and 


X.'a  are  ch.  roots  of  (7  A B Another  cheek  can  be  made  by  putting 
J M Ad  Ad 


P(pap)  * O(pxp).  Then 


y - X («p) 


Is  e central  chi-square  variable  with  wp  d.f.  and  G(0)  « G^(0).  Hence  we  get 


T - C((pts»fl)/(«|H-2))X^  + (p-nfl)X^]  + 0(n’^) 


2 2 

where  X > x (ap) . This  agrees  with  Ito's  result  [ 5]  to  the  Indicated  order. 


4.  AN  ASYMPTOTIC  EXPANSION  FOR  THE  C.D.F 

’ — 1 tin  — 1 


OF  T « m Tr  S^S2  UP  TO  0(n  *) 

In  this  section  we  will  derive  an  asymptotic  expansion  for  the  c.d.f. 


-1, 


of  T to  0(n  ),  following  the  method  described  earlier.  Again  we  write 
PtU  Tt  Sj^*<  20)  - J Pr{«  Tr  < 20|S2)Pr{d^) 


e Pr{m  Tr  < 20). 


From  Eq  (7)  we  get 


Pr{a  Tr  < 20)  - G(0)  ♦ T + 0(n"^) 


rs  tu  St  ur 


G(0)  - f[h,(A"')3G'(0)  + 0(n‘^>. 
n 1 m 


(23) 


Let  F(20)  be  the  c.d.f.  of  T;  l.e.. 


F(20)  - Pr{m  Tr  SjS"^  < 20) . 


Upon  substituting  Eq  (19)  in  Eq  (23)  and  replacing  y by  28  we  get 


F(20)  - G(0)  - ^ i A, (20)^1  g^(0)  + 0(n'^. 
™-j.l  J P 


(24) 


Hence  we  have  the  following  theorem : 

11 


Thaor—  2.  Under  the  eesuaptlon  (1)  of  Theorem  1,  the  asynptotic  expansion 
fox  c d.£.  of  T la  given  by  Eq  (24). 


Special  caeea; 

1.  Upon  neglecting  all  texna  Involving  where  f^^  la  the  (IJ)  eleawnt 

of  F,  we  get  from  Eq  (24) 

F(20)  - G(0)  - 4 r A,  (20)0  g^<0)  + 0(n’^), 

J P 

where  A^'a,  J - 1,2,3  are  the  aaaie  as  In  Eq  (21),  and  this  agrees  with 
C-P  [ 1]  to  the  Indicated  order  after  soam  simplifications. 


2.  Again  we  put  F(pxp)  - O(pxp),  and  we  get 


F(2G)  - Gp(0)  - ^(p-mH)0  + 2((pfmH)/(m|>f2))0*]  g^(B)  + 0(n‘^), 
and  this  agrees  with  Ito's  result  fS  ] to  the  Indicated  order. 


5.  AN  ASYMPTOTIC  EXPANSION  FOR  PERCENTILES 
OF  T » m Tr  S^S'^  TO  0(n‘^) 

.2 

Here,  using  the  technique  stated  earlier,  we  obtain  the  terms  of  order  n . 
The  results  of  the  third  and  fourth  derivatives  of  G(6)  as  it  stands  are  not 
convenient  for  practical  use.  In  order  to  make  some  simplifications  we  assume 
that  terms  involving  f^jf^,  are  negligable,  where  f. . is  the  (i,j)  element  of 
the  deviation  matrix  F.  From  the  third  and  fourth  derivatives  given  in  [9] 
and  from  Eq  (12),  using  a technique  similar  to  that  of  [5],  we  get  after 
tedious  simplifications 

Bj  . b'  + 24(cj2>/2))  - J - 1,...,4.  (26) 

The  coefficients  cf^^  and  cj^^  are  given  in  Eq  (32),  but  the  b!'s  and  C.'s 
J J J J 

are  listed  below- 
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L B.y  g(fl)[G'(0)r' 

48n^  J«1  1 ^ 


8n 


where 


b{  - 7p*  + (-12i»fl2)p  + a«^-12BflX 

* (13p^+24p-ll*^47)/(Bp4-2), 

b'  - (4mp^+2(3«^+3»flO)p*+2(2«^43m^+17m 
+ 18)p  4-  4 (5m^+9»f2))/(mp4-2)^(itp+4), 

- 6(p-l)(p<-2)(m-l)(BH-2)/(«pf2)^(«lM)(ii|>f6X 

C.  » (F)  ((p-»fl)  ((ia^/8)p(p+l)  - (o^/8)p)  + (-(m/4)p. 
i 

<pfl)  (m^/4)p)(-(m/2)(p4l)  + (m^/2))), 

C,  - (F)((p-»fl)(-(«/8)(p+l)  + («/2)/p  + (3«‘/8)) 

+ (p4ofl)  ((m^/8)p(p+l)  - (■^/8)p)  + (-(a/4)p. 

(lH-1)  + (m^/4)p)(-(2a)/pV)/(iip+2)  - (a/2)  (-(m/2) 
(R4-1)  + (m^/2))X 

C,  - (F)((p-afl)(-(i>fl)/2p  • (B/8)(pfl)-B/p-(3a^/8))/ 
(apf2)  + (i>4iiH-l)  (- (b/8)  (pf  1)  + (b/2)/p -f  (3m^/8))/ 
(mpf2)  + (p-mfl)(-(B/8)(p+l)  + (a/2)/p  + (3b^/8)) 

+ (-b/2) (-(2m)/p  - (B^))/(apf2)  + (-(m/2)(p+l) 

4 (B^/2))((pfl)/4  + (m/4)p)/(Bpf2))^ 

C.  - (F)  ( (p-arf  1)  ( (p+1) /2p  + (b/8)  (pfl)  + (a/2)/p 
+ (B^/8))/(apf2)  (apf4)  + (p+mfl)  (-(pfl)/2p 

- (b/8)  (pfl)  - (b/p)  - (3a^/8))/(apf2)*  +(-m/2) 
(2(pfl)/p  + B((pfl)/Z  f (2/p))4«^)/(Bpf2)  (mpf4) 

+ ((pfl)/4  + (b/4)p)  (-(2b)/p  - m^)/(apf2)^) 
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tJn 


Cj  - (r)((i>Hirt-l)((p+l)/2p  + (a/SUpfl)  + (a/2)/p 
+ (m^/8))/(mp+2)^(Bpf4)  + <(pfl)/4  + (n/4)p) 

(2(p+l)/p  + m((pfl)/2  + 2/(p)  + (m^/2))/(iii>f2)^(mp+4). 

Now  we  subatituce  Eqs  (21)  and  (25)  into  Eq  (8)  to  get 

£i£2^“  ^ ■*■ 

A yJ  g (0)[G'(e)r^  f-l~  s bV. 

^ ^ 24n^  J-1  J 

gp(e)[G'(0>r^  " 2 C yJ  g (0)rc'(0)]"^  + 0(n"^)  (2) 

4n  j«l  J P 

where  A^'a  are  given  by  Eq  (20).  The  Bj  »s  and  C. 's  are  given  above,  and 
get  aa  a final  reaulc  the  following: 

IhS2ISB_3.  Let  nS^  and  n^  be  independently  distributed  W(m,p,B’^)  and 
-1 

W(n,p,A  ),  respectively,  and  let 

i ‘i  ■ •"'>  < 1.  1 ■ 1 p.  and 

(11)  terns  involving  be  negligable.  where  f^j  is  the  (1,J) 

element  of  F.  Then  the  aaymFtotic  expansion  for  the  percentile  of 
T - m Tr  ia  given  by  Eq  (28). 

6.  AN  ASYMPTOTIC  EXPANSION  FOR  THE  C.D.F. 

OF  T = m Tr  S^S"^  TO  ORDER  n"^ 

Here  we  will  derive  an  asymptotic  expansion  for  the  c.d.f.  of  T 
following  the  methods  described  in  the  previous  pages.  Again  we  write 

Pr{m  Tr  8^82^  1 26}  = Pr{m  Tr  S^S"^  ^ 20}  Pr{dS2) 

= e Pr{m  Tr  SjA  ^20). 
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From  Eq  (7)  we  have 


Pr{m  Tr  S s'^  < 29)  - G(0)  + - E o a 9 & 

n n tu  ft  ur 

G(0)  + ^ E GO)  + 

2n^  ^ fy^ft^ur^wx^yv®^®^  + 0(o  ). 


(29) 


Ut  F(20)  be  the  c.d.f.  of  T,  l.e.. 

P(20)  - PrU  Tr  < 20). 

Opon  flubatltutlog  Eq  (19)  and  results  from  [ 9]  we  get  after  tedious 

simplifications  the  following: 


F(20)  . G(0)  - 


4 

E 

J-1 


(20)^3  g^(6)  - 


4 

r E 

J-1 


Bj  <0)^3  8p(0)  + 0(n’^) 


and 


64C 


(1) 


+ 24T. 


(2) 
J •* 


(30) 


(31) 


where  the  coefflclenta  A^'s,  are  available  in  [9]  and  the  rest  of  the 
coefficients  are  listed  below: 


15 


" 2(3«  •3af4)p^  + 3(a^*2B^<f5a-4)p  - 8a^-fl2Bf4)  p 
3 2 

bj  - -2(3«p  +2(3b  +3i«-4)p^  -3(3»^-2m*-5iH4)p  - 8«^+12n+^) / (mpf 2)  9 

bj  - 4(3mp  - 2(3b  -3m>4)p^  - 3(3B^+2m^+ll«-4)p  - 40m^  - 36m  - 4)/ 
(mpf2)  (inpf4) , 

b^  - 24  (mp  -f2(m^-HBf4)p^  + (m^+2B^+21mf20)p  + (8B^+20mf 20)  > / 

(np+2)  vnpr4)  (mp4-6) , 

* (rX-(m/16)(p^4-3i>f4)  4 OaVSMpfl)  - m^/8), 

(1) 

^2  “ (£)X-3/2)(p‘-4-3p+4)/p  - (3m/2Hp)-0/p 

- (3ra^/4)(pfl)  + (3m'"/2)/p4.  (3m^/4))/(mpf2) , 

“ ‘'|^X(6m)(p+l)/p  f (3m^/2)(pfl)  + (6m^/p 
+ (3m^/2))/(mp4-2)  (apf4), 

^4  ^ ^ +3p+4)  /p  + (m/2)  (p^4-3pf4) 

(I8m)(pfl)  /p  4.  (3m^)(p+l)  + (6m*) /p 
+ («^))/(mpf2)(Bp+4)(tnp46), 

(2(2p^4>5pf5)/p  + 'n/2H2p^+5pf5) 

^ (8m)  (pfl)/p  4-  <B)<p^42p^4;5jH-2)/^  4. 

(p  +pf4)/p  4-  (m^/4)  (pSp44)  4-  (m^/2)), 

^2  ^ “ (£)((14)(2p*+5p+4)/p  4-  (2a)  (2p^4-5p4-5) 

4-  (60a)(pf])/p  4-  vftm;  (p>2p^4lp»2)/p  I 

+ (6m^)  (p+l)  4-  (n^l  (p^+2p^+)pi-2)  4-  (9*^) 

(p^+p44)/p  4-  <3b"/2)  (p^4i>f4)  + <3m'*))/(apf2)  , 

- -(F)((48)(2p*4-5p+5)/p  + (8m)(2p*4-5pf5) 

4-  (I92m)(pfl)/p  4-(28m)(p^4-2p^4-3pf2)/p 
4-  (24m^)(pH)  4-  (4m^)  (p^4-2p^4-3pf2)  4-  (24m^) 

(p*+p+4)/p  4-  (3B^(p^4trt4)  4.  (6m^)/(^2)  (mp44) 


(32) 
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1 


I 

i 


iiliiliillilillll 


and 


I.-, 

r 


i 


i 


I 

j 


= (F)((64)(2p^+5p+5)/p  + (8m)  (2p2+5p+5) 

+ (208m) (p+l)/p  + (32m) (p^+2p^+3p+2)/p 
+ (24m^)(p+l)  + (4m^)  (p^+2p^+3p+2)  + (20m^) 

(p^+p+4)/p  + (2m^) (p^+p+4)  + (4m^))/(mp+2) 

(mp+4) (mp+6) . 

Then,  we  have  the  following  final  form: 

Theorem  4.  Under  the  same  assumptions  of  Theorem  1,  the  asymptotic  expansion 

_2 

for  the  c.d.f.  of  T to  0(n  ) is  given  by  Eqs  (30),  (31)  and  (32). 

7.  NUMERICAL  STUDY  OF  POWERS  AND  ACCURACY  COMPARISONS 

To  show  hok'  the  accuracy  has  improved  by  introducing  the  terms  involving 

f..f,,, where  f.  . is  the  (i,i)  element  of  the  deviation  matrix  F as  well  as 
ij  kl  ij  ^ ■’  * 

_2 

terms  of  order  n , some  numerical  results  are  presented  in  this  section. 

Some  comparisons  may  be  made  from  Table  I between  the  exact  and  approximate 

powers  of  T when  ps3  and  m»4  for  n>=34  and  n>84.  Values  of  the  exact  powers 

in  the  table  are  taken  from  Pillai  and  Sudjana  [4].  Our  expansion^which  is 

given  above  by  Eqs  (30),  (31)  and  (32), is  used  for  the  computation  of 

-1  -2 

this  table  up  to  0(n  ) and  0(n  ).  To  illustrate  the  usefulness  of  the 

neglected  terms  in  the  C-P  approximation  [1]  values  using  their  approximation 

are  given  in  ( ) in  the  table.  It  may  be  observed  that  their  values  differ 

considerably  from  the  exact  while  our  improvement  leaves  practically  very 

little  error.  One  can  see  from  Table  I that  the  approximation  given  by 

Chattopadhyay-Pillai  to  order  1/n  is  not  very  good  even  after  adding 

terms  of  order  1/n'*.  Again  from  Table  I it  is  obvious  that  the  accuracy 

given  by  terms  of  order  n'^  is  not  enough  even  after  including  those 

_2 

^ij^kl  usefulness  of  terms  of  order  n is  also  considerable. 

Further  power  computations  have  been  carried  out  for  p>3  and  p*4  and 
presented  in  [9].  For  tabulations  of  powers,  the  upper  five  percent 
points  were  taken  from  Davis  [10]. 
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TABLE  I 


TOE  COMPARISON  BETWEEN  THE  EXACT  AND  APPROXIMATE  POWERS  OT  T-TEST 
FOR  p=3,  m=4,  a=0.05  AND  FOR  EQUAL  DEVIATION  PARAMETERS. 


Up  to  the 

<D 

order 

n*34 

n =84 

0.001 

0(1) 

0.013 

0.031 

0(n“^) 

0.036 

0.048 

CM 

• 

C 

O 

0.049 

0.0501 

Exact 

0.050 

0.0501 

0.150 

0(1) 

0.019 

0.043 

0(nh 

0.049  (0.052) 

0.063  (0.065) 

0(n'h 

0.064 

0.065 

Exact 

0.064 

0.066 

o.soo 

0(1) 

0.040 

0.079 

0(n’b 

0.091  (0.097) 

0.108  (0.115) 

0(n"*) 

0.107 

0.109 

Exact 

0.102 

0.109 

1.000 

0(1) 

0.087 

0.150 

0(n’b 

0.179  (0.178) 

0.188  (0.206) 

0(n’^) 

0.202 

0.189 

Exact 

O.J73 

0.189 

The  figure!  In  ( ) ere  coaputed  ueing  Chattopadbyay-Pillel  expansion. 
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SECTION  II 

ASYMPTOTIC  FORMULAE  FOR  THE 
PERCENTILE  AND  C.D.F.  OF  HOTELLING»S 
TRACE  UNDER  VIOIATIONS 


I 

I 


? 


? 


i 


1.  INTRODUCTION 

In  the  previous  section,  asymptotic  expansions  for  the  distribution 

and  percentile  of  the  statistic  T = m Tr  have  been  obtained  up  to 

2 

terms  of  the  order  1/n  , where  mSj^  and  nS2  are  independently  distributed 
central  Wishart  with  m degrees  of  freedom  and  covariance  matrix 
Ej,W(m,p,Zj) , and  with  n degrees  of  freedom  and  covariance  matrix  Z2, 
W(n,p,Z2)»  respectively.  Further,  denoting  the  non-centrality  by 
(F)  « TrF  « Tr(B  ^A-I),  where  = B and  Z^^  = A,  we  also  included 
terms  involving  where  f^^j  is  the  (i,j)-^  element  of  F,  which 

were  previously  neglected  by  Chattopadhyay  and  Pillai  [1].  In  this 
section  again  we  extend  the  work  of  Chattopadhyay  [11],  who  derived  an 
asymptotic  expansion  up  to  terms  of  order  1/n,  neglecting  terms 

for  c.d.f.  and  percentile  of  the  trace  statistic  when  mS^  has  non-central 
Wishart  distribution  with  ra  degrees  of  freedom,  covariance  matrix  Z^  and 
non-centrality  parameter  n,  W(m,p,Zj,n)  and  nS2  distributed  central  Wishart 
W(n,p,Z_).  The  extension  in  this  case  is  to  include  the  f.  .f.,  terms 
neglected  by  him.  It  may  be  noted  that  these  terms  were  found  to  improve  the 
expansion  in  the  previous  section.  The  results  are  helpful  for  the  study 
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of  the  violation  of  a)  the  asstnption  of  a coaaon  covariance  utrix  in 
the  MANOVA  test  based  on  the  trace  statistic  and  b)  the  normality  assump- 
tion in  testing  » £2*  For  ■ £^>  asymptotic  expansions  of  the 
non-central  c.d.f.  have  been  studied  by  several  authors  [12]  and  [13]. 

2.  THE  METHOD  OF  ASYMPTOTIC  EXPANSION 

The  notations  in  this  section  generally  follow  those  of  the  previous 
Section  and  other  papers  referred  to  earlier  [1],  [5],  but  additional 
notations  will  be  introduced  here.  The  method  herein  is  also  to  obtain 
an  asymptotic  expansion  for  the  percentile  of  T first,  and  use  it  to 
derive  an  expansion  for  the  c.d.f.  of  T,  where  T may  be  defined  as 
follows: 

Let  Z > be  a pxm  matrix  of  independently  distributed 

columns  vectors,tdiere  has  the  density  N(]j^,£^),  i « l,...,m.  Then 

e ® « 

we  may  define  T ■ Tr  Sl  ZZ  ■ T *]Sl  z,  where  nS,  is  distributed 

•2  *"•  -1-2  -i  -2 

independently  of  Z. 

Now,  if  is  replaced  by  B in  T,  then  Tr  BZZ*  is  distributed  as  a 

non-central  chi-square  with  mp  degrees  of  freedom  and  non-centrality 
2 

parameter  o>  , where 

- Tr  B»M'  - Tr  n 

•P  « « 

u ■ I*  0,  P ■ mp/2 

we  may  note  that 

Pr{Tr  BZZ'  < 9) 

— tm  — 

. .-“V2  j J.W  /'  .-x/2*, 

j-0  Ji2’^r(p*j)  0 
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vp 

2 

where  G (e,(i>  ) is  the  c.d.f.  of  non-central  chi-square,  with  ap 

®P 

2 

degrees  of  freedoa  and  the  non-centrality  paraawter  u . 

Let 


G(e)  - Pr{Tr  m*  ^ 6). 

As  a first  approxination,  for  large  n we  nay  replace  by  in  G(6), 
and  consider 

G(0)  - Pr{Tr  ZZ’  £ 0}. 

Furtheraore,  as  suggested  by  Ito  [ 5],  obtain  a function  h(§2)  of  the 
eleaents  of  S2  and  n large  enough  such  that 

G(0)  - Pr{Tr  ZZ*  < hty)  (33) 

and  then  write  h(S2)  as  a series  with  the  first  tera  being  a linear 
function  of  non-ceritral  chi-square  variables  and  terns  of  decreasing 
order  of  nagnitude. 

Now  Eq  (33)  can  be  written  such  that 

G(0)  • Eg^(Pr[Tr  S’^  ZZ»  < h(S2)/S2J}.  (34) 

By  using  Taylor's  expansion  it  is  possible  to  expand  Pr[Tr  ZZ'  ^ 
h(S2)|§2]  about  an  origin  p)»  where 

^ * (^»j)  " (!»•  ••»?)•  (35) 

Thus, 
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Pr{Tr  §2^  ZZ*  ^h(S2)|S2> 


{explTr(S2-A“^)a]Pr{Tr  ^Z*  < h(A'^)}, 


where 


l(pxp)  . j|^) 


and  £^j  is  the  Kronecker  delta.  Hence 


where 


G(e)  - e . Pr{Tr  AZZ*  < h(A"*)} 


8 - exp[-Tr  A‘^3]|I  - 

■ 1 ♦ ~ Te  e^  3 ..3  ♦ 0(n”^) 

n * rs  tu  St  ur  *' 


where  £ denotes  the  suBMtion  over  all  suffixes  r,s, each  of  which 


range  free  1 to  p.  Expanding  h(S2)  around  6 gives 
h(§2)  - e ♦ hj(§2)  ♦ h2(S2)+... 

where  h^(S2)  is  0(n*^).  Proa  Eqs  (38)  and  (39)  and  expanding 
h(§2)  around  h(A~^)  we  can  get 

G(0)  - [1  ♦ ^ K,Otu®st*ur  * 1 ♦ hj(A’^)D  ♦ 0(n’^)l 


Pr{Tr  AZZ*  < 0), 

m •• 


where  D ■ and  by  equating  teras  of  successive  order  we  get 

[hj(A‘^)D  ♦ i I®rs®tu*stW**^^’*’’^  ^Z'  < 0)  - 0.  (4 

For  the  purpose  of  evaluating  Pr{Tr  AZZ'  ^ 0)  we  will  use  the 

perturbation  technique  [ 8J . 


^ ! 


J ■ Pr{Tr(A'^  ♦ e)"^  ZZ*  < e> 


where  e(pxp)  is  a syHMtric  aatrix  sufficiently  close  to  O(pxp). 
Following  [11],  [ 5]  and  [13],  we  get 

J - |l  - Bxp[-wV2] 

Exp{(l/2)E  Tr(I  - XA)"^Q)G  (6,0) 

where  A - E-1,  E'G^(6,«^)  - G_^,  (6.«^)  and 
■p  mp^^T 

X - ♦ e)-^  - I 

_ /b”!,  ^b”1*x  /*''l»^ 


• (»■**-!)  - ^„(!‘*A)(*„A) 


idiere  A^*  is  the  pxp  eatrix  with  (i,J)-^  eleeent 

Also  by  Taylor's  theorea  J can  be  expressed  in  the  following  fom 

J ■ (1  ♦ Tc  6 ♦ »r  I®  ®*  ^ 3*  ♦...} 

® rs  rs  21  ^ rs  tu  rs  tu 

Pr(Tr  ^Z'  < d)  (44) 


Now, if  B"^A-I  ■ E such  that  |chj(F)|  < 1,  i > l,...,p,  upon  using  the 


notations 


Tr(A;J  A)  - (rs) 

Tt(A;J  A)(A;i  A)  - (rs|tu) 

Tr(P)(A;J  A)(a;J  A)  - (F|rs|tu) 

Tr(B"*A)(A;J  A)(B'^A)(A’^  A)  - (I+F|rs|l+F|tu) 
Tr(F^)  - (F^),  Tr(F*)  - (F^),...  etc. 
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•nd  substituting  X in  Eq  (43)  ten  by  ten  coqperison  between  the 
two  expensions  of  J,  Eqs  (44)  and  (43),  after  substituting  X will  give 
the  second  derivative  9.«.3.^Pr{AZZ'  < 0), which  can  be  written  in  the 
following  fon: 

i •)  - C4S) 

idiere 

aJ  - ^{-2(I+P|rs|tu)  ♦ (I+F|rs|l+P|tu) 

♦ 2(P|l4P|rs|tu)  - 2(P|l4-P|rs|l4P|tu) 

.2 

- 2(P|F|I+P|rs|tu))+  ^){(I^P|rs)(I+P|tu) 

* 2(F)(I«P|rs|tu)  - (p(I«P|rs|l«P|tu)  • 2(F)  (P|l^P|rs|tu) 

- 2(I*P|rs)(F|H|tu)  - (P^)(I*E|rsltu)) 

- (^)<(P)(H|t»)(H|tu)  ♦ (P)^(i*F|rs|tu))  . 

Other  Aj's  coefficients  are  available  in  [14]  ^ [15], 

3.  AN  ASYMPTOTIC  EXPANSION  FOR  THE 
PERCENTILE  OF  T « Tr  S"^  ZZ* 

Recalling  that  G_(6*m  ) is  the  c.d.f.  of  the  non-central  chi-square 

distribution  with  mp  degrees  of  fteedoa  and  non-centrality  paraneter 
2 

tt  ,we  nay  note  that 

Hence,  it  is  possible  to  rewrite  Eq  (45)  in  the  following  fom, 
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Again,  we  note 

9 e w gU 

Zo„(rs)  - p.  Zo^^Oy,(rs)(tu)  - p, 

“ " • (f)* 

“ " • ^(9)CP)  ♦ (OF)].. ..etc. 

As  a check  for  the  above  relationships,  let  F(pxp)  ■ I(pxp)  and 
O(pxp)  - I(pxp);  thus  y should  equal  to  p.  which  is  the  value  of 
^®st°ur^“^  * Siailarly  v and  w will  be  reduced  to  Ioj^Oyy(*“*|tu) 

equal  to  l/2p(p«l]:  With  the  help  of  the  above  relationships,  it  is  possible 
to  evaluate  the  AJ's,  J ■ 0.....6.  after  suning  over  all  subscripts. 

Jfg  Sg  tg  u« 

Now  by  using  Eq  (41)  and  the  above  coefficients  we  get 
-hi(A"^)D  Pr{Tr  AZZ'  < e) 

1 ^ 

* « jlo 

where 

•o  ■ >5>0>-P-i) 

aj  - -2«(»p-«2). 

aj  - ^(■♦p+1)-2(2«+p+1)«2  ♦ tr  0*. 

»3  ■ 2{(«+p+l)w^  - tr  5^), 

- tr  ^ 
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ERS.'  ffl! 


* f »# VfU 


3 * 0f***f6< 


The  above  can  be  sUplified  after  tedious  algebra,  and  it  can  be 
written  in  the  following  si^>le  expression. 


■ Jo 


«J  * 0,1, 


(47) 


Some  of  the  Ajj^'s  are  listed  below  and  the  rest  are  listed  in  [14] 


Aqi  - -(1/4)[2[(P)^+(f’)1(P*1)4[(P)^(f2)*2(f2)(P)+2(P*)]), 

Aq2  - (1/8)[((P)p-2(P)2-S(E)(P^)-(P^)^1(P*1) 

- ((P)’>(F)(P^)+6(P^)44(F*)33, 

/•OS  • -(1/8)[(1/2)[(P)V(P)^HP*1)*(P)P*2(F)^+(P)(P^)), 

^10  • <-5(9)l(P^)*(P*)l  - 5l(0P^)  ♦ 

Aji  ■ {(24(F2)4l8(p5)*(Q)(p5)  ♦ [(3/2)p*S(F*)*3(P*)3.  (OP) 

♦ (-(3/2)p  ♦ 3(P)*3(f2)3(OF^)  ♦ [3p*3(P)3(OP’)](p*l) 

♦ [(9/2)(f2)^(4S/2)(F)(F^)*19(F)(p’)3(0)  ♦ [(3/2)(F)^ 

* ^ t mr  mt  w « 

♦ (9/2)  (P^)*3(F)  (F^)*4(P’)3  (QP)*[6*21(F)+(3/2)  (P)^ 

♦ (9/2)  (F^)3  (OF^)*(18(F)*543  (OP’)  *36(09^) +12  (F)^ 

• mm  m mm  mm  m 

♦ 12(P^)^18(F)(F^)+18(F^))/(12), 

m mm  m 

A,.  - <[-(P)P*3(FV8(P)*+9(F)CF^)*[(F)(p2)-(p2) 

JLm  * * m mm  m m m 


26 


i 


♦ (P)^l  (?)/2*[-(P)lH2(P)^S(P)(P^)*(p2)l  (0F)/2 
♦I(P)lH(P)*l(eF^)l.  Cp*l)*[-2(P)*(P^)-(l/2). 

(F)  (P^)+(P’)*2(p2)  (P)^^(P)*]  (9)*I-p*2(F) 

♦ (9/2)(P^)*2(P)^+(P)(P^)+2(^)+(F)V2]  (OF) 

♦ Ip*12(F)+(3/2)  (F)^+2(F^)1  {QP^)+6(QP*)  (F) 

* « « «»«»« 

♦ 4(P)420(F^)+12(F^))/(8). 

Ajj  - {[9(P)V(F^)(F)^*(l/2)(n)(F)^t(F)^(F)^l.  (3(QF)/2)]  . (p+1) 

♦ [3(F)  ^+(1/2)  (F)^l  (Q)+(3(F)p*12(P)243(P)  (F^) 

♦ (1/2)  (F)*)  (nP)*3(QF^)  (P)^+18(F)p*36(P)^ 

♦ 18(F^)(P))/(48), 

Ajo  ■ <-(9r*)(p*l)^[4(P)442(p2)*36(F’)Hn) 

♦ (4-3(0)  (F)-3(0)  (p2)J  (0F)*(38-(0)  (F)]  (OF^) 

- 3 (gp)  *+28  (OF^)  -4  (OF^)  ♦ (QFO* ) ♦ (0 ' TO) 

♦ (n'F^n)-4(0F)(QF^))/(4). 

Aji  - {[-8(F)-40(F^)-24(f’)*(0)[4(F)-2(F^)-(16/3)(f’)] 

♦ (-2p-20(F^)-16(F^)l(0P)*(4p-20(P)-16(P^)](0P^)-(16) 

IP+(P)]  (QP’)^(-P/2*(F)*(F^)]  (0P)^[2p+2(F))  (OF)  (OF^) 


* (?)  (??)  (?^)1  (P+1)-(«(F)^*1S4(P)  (P^)*^P)  (F*)+24(p2)^]  (8) 

- [8(F) -16*10 (F)  ^*36 (F^)+|i{P*) ♦16(F)  (F^)J  (OF) 

- [64*98(P)*S2(F^)*8(P)^) (OP^)- (312*96 (F)] (OF*) 

- 192(flF^)*(4*6(F)*(3/2)  (P^)*(l/2)  (F)^]  (OT)^ 

+ [-4*6 (F) ^*7  (F)  (F^)]  (0)  (OF)  -20 (F)  ^-20 (P^)  -24  (F)  (F^) 

* 24 (OF') -24 (F*)* [-2(P)*(F)^1  (OTn*)*[6(p 

* (F^)  ] (n • FO) *[14  (F) *(F^)  J (R*  F^O) * [28*6  (P)  ) (OF) . (QF^) 

* 12(nP)  (nF^)*[-2(F)*(P^)]  (n'Q)+8(F)  (fl'P^fl) 

* 4(n)(flF^)*4(nF^)^)/(16). 

A22  - ([-4(F)p-40(F)2-36(F)(f2)-12(p2)P-[8(F)(P^) 

* 2(P)^1  (fl)*[4(F)p-20(F)^-24(F)  (F^)-8(P^)p]  (OF) 

- 16[  (F)p*(F)2]  (W^)*[(F)p*(F)^1  (9F)^(9)  (OF)  (F)^] . (p*l) 

* [24  (F)  - 16 (F^)  -IS  (F) ’-32 (F*)  (F) (n) * [ 1 2p-8 (F) ’-32  (p’) 

- 24(P)-84(F^)-26(F)^-16(F)(F^)]  (OF)-[16p*200(F) 

^ ••  ^ mm 

* 32(F^)*24(F)^1  (nF^)*[p*10(P)*(F^)*(P)^]  (OF)^ 

• ^ ^ mt  m 

- 32(F)-88(F^)-12(F)’-12(F)(F^)-48(F’)-32(F)(nF')  j 

- 64(F)  (flF’)*(Q'0)  (F)^*(nPn*)  (F)^*(0'FQ)  (F)^ 

m m m ^ ^ mt  « w «»  iw 
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♦ (? ' F^O)  (F)  ^+4  (OF)  (P)  +4  (QF)  (fiP^)  (F)  )/  (32) , 

Ajj  ■ -n[(3/2)(F)V(S/2)(P)’)(flP)*a/2)(0)(F)’].(p4l) 

♦ [S(P)p*12(P)^+3(F^)  (F)+(l/2)  (P)’l  (QF)  + [3(F)  ^+(1/2) 

(F)^](9)*3(P)^(?F^)}/(12).  Mid  the  rest  available  in  [76]. 

As  an  iaaediate  result  h(§2)  can  be  expressed  in  the  following  Banner 

h(Sj)  . e-[]j 

* n AjG^2j(®»‘^>H6'(e)]"^  0(n"^).  (44) 

Recall  that  6 is  the  appropriate  percentile  of  the  linear  function  of 

a non*central  chi-square  variable  of  the  fora  Y ■ X X4X. (a»ii)  )»  the 

j-1  J J 

X.  *s  are  the  characteristic  roots  of  AB  and  G(0)  is  the  c.d.f.  of 

j 

Y in  terns  of  the  percentile.  Finally,  we  can  state  the  following 
theorea: 

Theorea  5.  Let  Z • (*jf***|,)  ^ » pxa  randoa  aatrix  of  independently 
distributed  coluan  vectors,where  Zj  has  the  density  N(yj,  ^j»B"^),  and 
n§2  distributed  central  Wishart  W(n,p,^2  ■ Under  the  assuaption 

that  B ^A  ■ I+F  and  |Chj  (P)|  < I,  i ■ l,...,p,  an  asymptotic 
expansion  for  the  percentile  of  T is  given  by  Eq  (48). 

The  following  are  special  cases  of  Bq  (48). 

Case_l.  When  terns  involving  ere  negligible,  where  f^^  is  the 

(i,J)  eleaent  of  F,  tens  like  (F)*,  (F*)  and  (n)(F*)  can  be  dropped. 
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Consequently disappears  and  up  to  A^are  reduced  drastically 

and,  finally,  Eq  (48)  agrees  with  the  result  of  Chattopadhyay  [11]  to  the 
indicated  order. 

Case  2.  Under  the  equality  of  the  two  covariance  aatrices,  the  deviation 
natrix  is  zero.  Putting  (F)  ■ 0 in  Eq.  (48),  we  get  Eq  (6.4)  of 
Siotani  [13]. 

4.  AN  ASYMPTOTIC  EXPANSION  FOR 
THE  C.D.F.  OF  T - Tr  S"^  ZZ* 

In  the  following  an  asy^>totic  expansion  for  the  c.d.f.  of  T to 
0(n"^)  is  derived  by  using  the  method  described  earlier.  Also, 

it  is  possible  to  write 

Pr{Tr  Sl^ZZ'  < 6)  - / Pr{Tr  Sl^ZZ'  < els,)Pr{dS,) 

-2  -2  --  — '-2  -2 


- 8 PriTr  WZ*  ^ 0} 


where  6 is  given  by  Eq  (39) . It  follows  that 

Pr(Tr  < e)  - G(0)  - ^hj(A"^]G' (0)  ♦ 0(n’^) 

Under  the  assumptions  of  Theorea  5 , we  get 

1 1 ^ 0 
Pr<T»  Sj*22'  < »)  . cm  * fe  I > 


where  the  a. 's  and  A.'s  are  presented  earlier  and  G(6)  and  G (e.u 
J J mp 

are  defined  earlier.  Then  it  follows: 
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Theorem  6.  Under  the  assmupt  i.  ;is  of  the  previous  theorem,  an  asymptotic 
expansion  for  the  c.d.f.  of  T is  given  by  Eq  (49). 

Similarly  we  can  get  the  two  special  cases,  as  we  pointed  out  in 
the  previous  pages. 

5.  NUMERICAL  RESULTS 

The  expansion  given  by  Eq  (49)  has  been  used  here  to  compute 
the  powers  of  the  test  when  the  departure  from  the  null  hypothesis 
occurs.  The  following  table  shows  these  powers.  For  this  tabulation, 
the  upper  five  percent  points  were  taken  from  Pillai  and  Jayachandran 
[16]. 
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TABLE  II 


POWERS  OF  T TEST  UNDER  VIOLATIONS  FOR  p«2,  m»3,  o=0.05 
WHEN  TOE  DEVIATION  MATRIX  HAS  EQUAL  DEVIATION  PARAMETERS 


“l 

“2 

(F) 

Up  to  the 
order 

n * 83 

.00001 

0(1) 

.0379288 

0(n"^) 

.049277  (.049273) 

0 

.00001 

.00015 

0(1) 

.037942 

0(n‘b 

.049296  (.049344) 

.005 

0(1) 

.038399 

0(n‘b 

.049433  (.04977) 

.00001 

0(1) 

.03793 

0(n‘b 

.049279  (.049276) 

0 

.0001 

.00015 

0(1) 

.037945 

0(n“b 

.049300  (.049347) 

.005 

0(1) 

.038403 

0(n'b 

.049636  (.05182) 

.00001 

0(1) 

.038085 

0(n‘b 

.04944  (.049455) 

0 

.005 

.00015 

0(1) 

.038098 

0(n"b 

.049437  (.049475) 

.005 

0(1) 

.038557 

0(n'b 

.0499188  (.05200) 

The  figures  in  ( ) are  conputed  using  Chattopedhyay  expansion  [11]. 
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